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Abstract 

The I-V characteristics of a Josephson junction shunted by an Ohmic resistor 

shows sharp peaks when levels in neighbouring wells are crossing. We consider 

the shape and size of these peaks using a double-well model where the wells 

are given in parabolic approximation. The friction of arbitrary strength is 

included by the help of a bath of infinitely many degrees of freedom. 
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I. INTRODUCTION 

The past decade has seen a considerable interest and remarkable activity in an area which 
presently is often referred to as macroscopic quantum mechanics. Specifically, one has been 
interested in quantum phenomena of macroscopic objects Jl|. 

In particular, macroscopic quantum tunneling (quantum decay of a metastable state), 
and quantum coherence || have been studied. Soon, it became clear that dissipation has 
a profound influence on these quantum phenomena. Phenomenologically, dissipation is the 
consequence of an interaction of the object with an environment which can be thought of as 
consisting of infinitely many degrees of freedom. Specifically, the environmental degrees of 
freedom may be chosen to be harmonic oscillators such that we may consider the dissipation 
as a process where excitations, that are phonons, are emitted and absorbed. Such a model 
(Caldeira-Leggett model) has been the basis in Ref. || where the influence of dissipation on 
tunneling has been explored. 

As far as quantum coherence is concerned, the most simple system is an object with 
two different quantum states: it is thought to represent the limiting case of an object in a 
double-well potential where only the lowest energy states in each of the two wells is relevant 
and where the tunneling through the separating barrier allows for transitions that probe the 
coherence. Since a two-state system is equivalent to a spin-one-half problem, this standard 
system is often referred to by this name. In particular, with the standard coupling to a 
dissipative environment made of harmonic oscillators, it is called the spin-boson problem 
which has been studied repeatedly in the past 0,|J. 

Level quantization and resonant tunneling have been observed recently in a double- 
well quantum-dot system. However, the influence of dissipation was not considered in this 
experiment. On the other hand, it seems that Josephson junctions are also suitable systems 
for obtaining experimental evidence pertaining to macroscopic quantum effects. In this 
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context, evidence for level quantization ||, [|7| and for quantum decay || have been obtained; 
however, macroscopic quantum coherence has evaded observation so far. 

Specifically, a Josephson junction may be characterized by a current-phase relation 

/(</?) = I j sin <p , (1) 

where the phase <p is related to the voltage difference U by 

hp = 2eU . (2) 

Therefore, the phase of a Josephson junction shunted by a capacitance C and biased by an 
external current I x obeys a classical type of equation of motion 



with the mass 



and the potential energy 



M^gc (4) 



V((f) = [I J cos (p + I x cp) . (5) 



A widely discussed model of a dissipative object is the one where the Josephson junction 
is also shunted by an Ohmic resistor R. In this case, the classical equation of motion @ 
has to be replaced by 

M<p= j^-Tlfp, 

dip 

The model of a dissipative environment according to the above specification has been dis- 
cussed by Caldeira and Leggett 0. 

The potential energy V((p) of Eq. displays wells at tp ~ 2nn with depth shifted by 
an amount A ~ (2irK/2e)I x . If the wells are sufficiently deep, one needs to concentrate only 



on transitions between pairs of adjacent wells. Thus, one arrives at the double well problem 
mentioned above. 

The analysis in this paper goes beyond the limiting situation where only the lowest level 
in each of the two wells is of importance. Roughly, this is realized when the level separation 
h(2Ej/M) 1 ' 2 ~ ftehlj/C) 1 / 2 is smaller than or comparable with A. In particular, we will 
concentrate on resonance phenomena which are expected to show up whenever two levels in 
the adjacent wells happen to cross when the bias current I x , that is A, is varied. 

For such values of the bias current, there appear sharp asymmetric peaks in the current- 
voltage characteristic of the Josephson junction. This phenomenon has been studied by one 
of us (A.S., together with A.I. Larkin and Yu. N. Ovchinnikov) || within the standard 
model in the one-phonon approximation. For bias currents that correspond to crossings of 
the next and next nearest levels (e.g. ground state in the left well and the first or second 
excited state at the right side), it is possible to neglect processes in the reverse direction 
provided that the temperature is sufficiently low. Thus, the restriction to a double well 
system receives additional support. 

The transfer of the object from the left to the right potential well is accompanied by 
the emission of an infinite number of phonons. Therefore, in a recent paper by Ovchinnikov 
and one of us (A.S.) [10|] the fact is taken into account that in the resonance region, the 



contribution of phonons of small energy is important as well as the contribution of resonance 
phonons with energy equal to the distance between levels in the wells. 

The above analysis relies on a perturbation theory for weak dissipation. Moreover, the 
calculations in that article are rather intricate. Therefore, we want to reconsider the problem 
in a more transparent way and thereby get rid of the restriction of small dissipation. 

II. DESCRIPTION OF THE MODEL HAMILTON! AN AND ITS EIGENSTATES 

The model consists of a particle, called "'object'" (coordinate R\), which is coupled (in 
the sense of Caldeira and Leggett 0) to a "'bath'" of harmonic oscillators (coordinates Rj). 
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We shall use the conventions j G {2, . . . , N} for the bath oscillators and k G {1, . . . , N} for 
the indices of all coordinates in the model. The double-well potential is approximated by 
two parabolas about the minima of the two wells as it is depicted in fig. [3]. 

The phase (p of the Josephson contact then corresponds to the object coordinate R\ of 
the model, and the voltage U is related to the tunneling rate J by 2eU = <Jp = 2-kJ. As it 
has already been remarked, the current I x is proportional to the bias A of the two wells. 
Thus, calculating the transition rate for different values of the bias A is equivalent to the 
determination of the I-V characteristics. 

Specifically, we want to write the Hamiltonian of the model in the form 

k j 

^l^E^iafi^. (7) 

The states for the two situations "'object in the left well"' and "'object in the right well"' 
will be denoted by \A L , L) and |A_r, R), respectively. If one projects onto the eigenstates |n) 
of the one- dimensional harmonic oscillator and takes into account the shift of the wells, one 
arrives at the following decomposition (ip n (R) = (R\n)): 

(n L ,{Rj}\A L ,L) =: J dRupn^ + aifattRk}), (8) 
(riR^RjjlA^R) =: J dR x ^„ H (fli - a) <f>% R {{R k }) . (9) 

The situations "'object on the left"' and "'object on the right"' differ only by the shift 
and the bias of the wells. Therefore, one can find a unified representation by noting that 
(J)^({R k }) = & A ({R k + a}) and <j>it({Rk}) = ^\({Rk - a}). The eigenstates <3> A are defined 



by the relations 

$ A ({i4}):=({i4}|A>, (10) 

H \A)=E A \A), (11) 

ffo = ^ + ?E - ^i) 2 + y ^R\ . (12) 

k j 
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Thus, it follows from eqs. @ and respectively, that 

(n L , {RjjlAL, L) = (n L ,{R j }\exp(iaJ2Pj)\^L>, (13) 

i 

(n R , {Rj}\A R , R) = (n R , {Rj}\ exp (-iaJ^Pj) \ A R > , ( 14 ) 

3 

where we have used the shift property of the momentum operator p. 

The coupling of the two wells is taken into account by means of a tunneling Hamiltonian 
Ht which we represent in the form 

(A l ,L\H t \Ar,R) = f d{R,} ]T T^iA^Lln^iR^inn^R^lA^R). (15) 

n L n R 

Using again the momentum operator, one can write 

\x) (x'\ = e ^ x '- x) \x') (x'\ = e if,{x '- x) 6(x' - x) . (16) 
From this, we conclude that 

(A L , L\H T \A R , R) = £ T nLnR f dR x dR' 1 d ^-^ nR {R! 1 )< PnL {R x ) 

n L n R J Z7r 

x (A L , L\e iMR i- Rl) e iQ{R i- Al) \A R , R) . (17) 



III. TRANSITION RATE 

The net transition rate from the left well to the right one is then in second order pertur- 
bation theory given by 

J = 2nZo 1 J2 \(Al,L\H t \A r ,R)\ 2 5(E Al -E Ar + A) 

x [e-^ - e^ A «] , (18) 

where Z = Trexp (—(3H ). The 5-function may be written in Fourier representation, and 
the fact that the Ea are eigenenergies of Hq serves us to incorporate the energy conservation 
into Heisenberg time-dependent operators A(t) = exp (iH t)Aexp (—iH t), i.e. 
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(A L , L\H T \A R , R) 5(E Al - E Ar + A) = 

= / dte iAt (A L \e- ia ^ pk{t) H T (t)e~ m ^P k{t) \A R ) . (19) 

Then, collecting our results from eqs. (0), (0), and (|19"D, we arrive at the expres- 

sion 



X 



Tr {e _/3 ^° e - 2ia £fc&(*) e ipi(*)(-Ri--Ri+2a) e -iQ(«i(t)--Ri) 



x e + 2ia J2kPk giPl (fix-/?! -2a) g-jQ^!-^)! ^ ^20) 

Let us now use the relation 



e 



-i(H a +W)t = e -iH t f e -i J* dt> W{t') ^21) 



which holds for t > when T is the time-ordering operator and for t < when the anti 
time-ordering is used. If we define (A) := Tr exp (— f3H )A/Z , we can write the following 
result for the transition rate: 



X 



dQdQ 



(2n 



dRdR'dRdR' e*Q^+*Q^ 



x ip nL (R) <p* R (K - 2a) ip± L {R ) ^ R (R - 2a) 



x (f exp [-iQi?i(t) + ipi{t)(R' -R) + 2imQ 2 a J* dt' R^t') 

-iQR^Q) + ipi{0)(R' -R)}). (22) 



We are now in the position to make use of the fact that the Hamiltonian is quadratic in 
all coordinates so that we can evaluate exactly 

^f e i J df „(f)fli(f)) = e ~UI dt'dt" v (t')D(t',t")v(t") ) (23) 
D(t',t") := —i(fR 1 {t')R 1 (t")) . (24) 



By comparison with the last two lines in eq. (|22|), the function i](t') is given by 

v (t') = -Q5(t' -t)- Q5(t') + 2mQ 2 a[Q(t') - G(t' - t)} 

+m(R - R')8'{t' -t)+ m(R - Tt)8'(f) . (25) 

is meant to represent the step function. The derivatives of the 5-function arise from a 
partial integration of terms containing p(t) = mdx(t)/dt. Note, that these act only on the 
coordinates but not on the step functions which arise due to the time ordering. 

Moreover, the degrees of freedom of the bath can be integrated out in the usual way |2| 
leading to a dissipative influence on the object. One is then lead to the following form of 
the Fourier transform of D(t, t') = D(t — t')\ 

^ x D R (u) D R (-co) 

1 — exp (— Tiuo/k B l ) 1 — exp {Tiuj/k B l ) 

(D R y\u) = m[{uj + iO) 2 - Q 2 } + ir)UJ , (27) 

where we will use a spectral density J{uj) = tjuj, < ui < uo c for the bath oscillators (see also 
[ iT| and appendix |A|). 

From eqs. (B2I) and (EBf) one can conclude that the integrations with respect to 



Q, Q, R, R', R, R can be done exactly as only Gaussian integrals are involved (note, that 
the eigenstates of the harmonic oscillator are Gaussian functions and derivatives of these, 
respectively). Therefore, for given riL,n R ,n L ,n R , one has to perform a 6-dimensional Gaus- 
sian integral. The only difficulty lies in the t-integration, which we have done numerically. 
Note furthermore, that without dissipation (i.e. rj = 0) the tunneling current J is indepen- 
dent of the shift a, as one can show explicitly. 



IV. TRANSITION RATE AT RESONANCE 

Specifically, let us restrict ourselves to the case of zero temperature where only the lowest 
level of the left well is expected to contribute to the transition rate. Moreover, we will assume 
that the bias A « Q and that the dissipation leads only to a moderate broadening of the 
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energy levels of the oscillator. Then, one needs to consider only the levels til = n-L — and 
n R = n R = I. 

The calculation of the transition rate for this case is deferred to appendix [FJ. We want to 
mention, however, that the cutoff frequency uj c of the bath spectrum is important to avoid a 
divergency of D(t) for t — > 0. In fig. ^|, we show the results for the transition rate calculated 
for different values of the dissipation 77 and the shift a. There, we have used dimensionless 
variables which means that the frequencies and lengths were normalized to their character- 



istic scales fl = \j2eIj/TiC and = ^jh/MQ, respectively. The dimensionless parameters 
are then given by a* = na for the shift and rj* = 1/(RCQ) for the viscosity. 

In a typical experimental situation, one has Ij ~ InA for the Josephson current and 
C ~ IfF for the capacity. Thus, we have k _1 pa 4.2, which means that a* ~ ttk pa 0.75. 
The dimensionless viscosity is given by 77* ~ 18kQ/R. Note, however, that the transition 
rate to be observed experimentally will be very small due to the exponentially small factor 
|T 01 | 2 /(M}) 2 §. 

For small values of the dissipation or rather of a = 2i]a 2 /7r, the perturbation theory [T^ 
yields the expression 

J = 2|T 01 | 2 Re (e™ a (^) 2a -^-), (28) 

where 7 = 77/ (2m) and e = A — Q. When normalized to the same value of the peak 
maximum, the curves calculated from the above formula (|28|) are nearly identical with those 
we obtained from our theory (cf. fig. |3|). The different normalization stems from [0 the 
way the divergencies were removed in ||10|| . 

V. CONCLUSION 

In summary, we have studied the transition rate in a double-well system under the 
influence of a dissipative environment. The rate shows asymmetric peaks when the levels 
of the wells approach each other. Thus, we have calculated the I-V characteristics of a 



9 



Josephson junction shunted by an Ohmic resistor for resonance conditions. We have modeled 
the system by two parabolic potentials coupled to a bath of harmonic oscillators. The 
coupling of the wells has been realized by a tunneling Hamiltonian. Moreover, the analysis 
has taken into account that the transfer of the object between the wells is a multiphonon 
process. In addition, our calculation has the advantage that the bath coordinates do not 
appear explicitly any more when the thermal average is taken. 
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APPENDIX A: GREEN'S FUNCTIONS FOR THE HARMONIC OSCILLATOR 

IN A DISSIPATIVE ENVIRONMENT 

From the equations of motion of the operators Rkif) and pk{t) in the Heisenberg repre- 
sentation, one can immediately infer the form of the retarded Green's function D R (t — t') = 
-ie(t-t')([x(t),x{t')}): 

(D R )-\u) = m[(u + *0) 2 - fi 2 ] - - / dco'— ^ T2 ^ (Al) 

7T JO \0J + 1\}) Z — LO 1 CO 1 

= m[(uj + iO) 2 - VI 2 } + ir/uj (A2) 

if the spectral density 0] is given by 

J{uS) = — mUj 5{uj — ojj) = Tju , < uj < uj c . (A3) 
2 j 

In order to avoid divergencies of D(t — > 0), it is necessary to introduce explicitly a cutoff 
for J(uj) with the characteristic frequency uo c . We used a Lorentzian cutoff which permits 
to evaluate the integral in (|A1|) exactly. For i]/(2m) C H C w c we obtain 

mD{t) = i Im (Ae~**) -ilmA e-^-vM 
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1 



A(g(-nt)+g{nt)) 



Im 

7T 
2 

+ — Im A Re g (iu c t — irjt / m] 



7T 



where we have introduced 

Q + ioj r 

A 



2Vt (ft + iu c — ir]/m) 
VI = Tt-ir)/{2m) , 



vl = n ^1 + 
rj — r] I 1 — 



n 2 \ 



The function g(oo) is related to the exponential integral 

Finally, the correlation function D (t — t') = —i([ct(t) , x(t')}+) is given by 

^ = 2.coth^Im^, 

so that we have 

2D(t) = D R (t) + D R (-t) + D K (t) , 
which immediately leads to eq. (|26|) . 

APPENDIX B: EVALUATION OF THE 6-DIMENSIONAL GAUSSIAN 

INTEGRAL 



First, let us write the transition rate ( p2|) in the form 

J = 2 Re / dte iAt K(t). 
Jo 

Also, we introduce new integration variables x = (v, q, u, v', q', u') such that 

R = (u - v - u' + v')/4 + a , 
R' = (u + v - u' - v')/4 + a , 
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R = (y + v + u' + v')/A + a , 
r' = (y, - v + u ' - v')/A + a , 
Q = (q-q')/2, 
Q = (q + q')/2. 

In addition, the eigenstates of the harmonic oscillator are given by 

V , ( J R) = v ^7T- 1 / 4 e -^ 2 / 2 , 

(px{R) = \f2Rnip Q (R) , with 
k = VmQ . 



Thus, the expression for K(t), t > 0, reads 

dqdq' 



K(t) 



x 



4tt 2 
d 2 



dudu'dvdv' \ T m \ 



16n 



dada' 



- exp (— xMx + bx + c) 



a=0=a' 

The only nonzero elements of the symmetric 6x6 matrix M are given by 



2 • 2 

k vm 
M n = ~< 7T 



(D(O)-D(t)), 



M 22 = -(D(0) + D(t)), 



K 



M 33 = — = M 66 



im 



Mia = M 21 = ~^D(t) = -M 45 = -M 54 , 
M 23 = M 32 = -| = M 56 = M 65 , 



M 44 = 



2 • 2 

k vm 



(D(0) + D(t)), 



M X = ~{D(0) - D(t)) . 

Thus, it can be decomposed into two 3x3 blocks. Furthermore, the vector b reads 

b =^(1,0, 1,1, 0,1) -^(1,0, 1,-1, 0,-1) 

+ ( K K 2 a-2im 2 Vt 2 a[D{t) - D{0)], ia + 2imQ 2 a J* dt'D(t'), 0,0, 0,o) , 
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and 



c = —2n 2 a 2 — ina(a — a') — i2m 2 VL A a 2 / dt\dt% D(t\ — t 2 ) — i2mVL 2 a 2 t . 

Jo 



(B9) 



The Gaussian integral in eq. (|B6|) can be evaluated exactly and eventually, we obtain 



K(t) = W(t) e 



s(t) 



K 



W{t) = — |T i| 2 (detM) 



-1/2 



, db , _ , i . db 



, , . , i . db dc 



S(t) = {hiM-^b + c} 

I ^ J a=0=a' 



a=0=a' 



(BIO) 

(Bll) 
(B12) 



In the case of zero dissipation, this expression is indeed independent of the shift a, as it can 
be shown exactly. 
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FIG. 2. Transition rate J*/!^] 2 as a function of the bias A* for T = (dimensionless, lo* = 60) 
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FIG. 3. Comparison of J*(A*)/|T * 1 | 2 (solid line) with perturbation theory (crosses) (dimen- 
sionless, rj* = 0.01, a* = 4, u* = 80, T = 0) 



17 



